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Abstract 

Introduction: Fixed orthodontic appliances, despite years of research and 
development, still raise a lot of controversy because of its potentially destructive 
influence on enamel. Therefore, it is necessary to quantitatively assess the condition 
and therein the thickness of tooth enamel in order to select the appropriate 
orthodontic bonding and debonding methodology as well as to assess the quality 
of enamel after treatment and clean-up procedure in order to choose the most 
advantageous course of treatment. One of the assessment methods is computed 
tomography where the measurement of enamel thickness and the 3D reconstruction 
of image sequences can be performed fully automatically. 

Material and method: OCT images of 180 teeth were obtained from the Topcon 
3D OCT-2000 camera. The images were obtained in vitro by performing sequentially 
7 stages of treatment on all the teeth: before any interference into enamel, polishing 
with orthodontic paste, etching and application of a bonding system, orthodontic 
bracket bonding, orthodontic bracket removal, cleaning off adhesive residue. A dedicated 
method for the analysis and processing of images involving median filtering, 
mathematical morphology, binarization, polynomial approximation and the active 
contour method has been proposed. 

Results: The obtained results enable automatic measurement of tooth enamel 
thickness in 5 seconds using the Core i5 CPU M460 @ 2.5GHz 4GB RAM. For one 
patient, the proposed method of analysis confirms enamel thickness loss of 80 (jm 
(from 730 ± 165 pm to 650 ± 129 |jm) after polishing with paste, enamel thickness loss 
of 435 pm (from 730 ± 1 65 liiti to 295 ± 55 Lim) after etching and bonding resin 
application, growth of a layer having a thickness of 265 liiti (from 295 ± 55 liiti to 
560 ± 98 pm after etching) which is the adhesive system. After removing an 
orthodontic bracket, the adhesive residue was 105 liiti and after cleaning it off, the 
enamel thickness was 605 [im. The enamel thickness before and after the whole 
treatment decreased by about 125 |jm. 

Conclusions: This paper presents an automatic quantitative method for the 
assessment of tooth enamel thickness. This method has proven to be an effective 
diagnostic tool that allows evaluation of the surface and cross section of tooth enamel 
after orthodontic treatment with fixed thin-arched braces and proper selection of the 
methodology and course of treatment. 

Keywords: Image processing, Orthodontic, Measurement Automation, Segmentation, 
Tooth enamel 
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Introduction 

Enamel is a tissue covering the dentine within the tooth crown. Its thickness in specific 
areas may exceed 1400 um. Enamel is composed of about 98% of inorganic matter in the 
form of di-hydroxyapatite crystals. This chemical composition gives it unusual hardness. 
Organic compounds and water constitute the remaining 2% of enamel [1,2]. Orthodontic 
treatment concerning a significant part of the population is related to enamel when it comes 
to installing fixed braces. In this case, there occur changes in the structure of enamel after 
treatment with fixed orthodontic appliances which depend on the used orthodontic hooks. 
There are different methods of attaching fixed braces and treating the tissue after hook re- 
moval. There is also different efficiency of cleaning a tooth surface off the composite mater- 
ial residue after the removal of orthodontic hooks when using selected methods of enamel 
treatment [3]. 

For example, there are techniques of adhesive attachment of orthodontic hooks to the en- 
amel surface which enable rapid development of orthodontic treatment with fixed braces 
[4-6]. This method is continuously being developed and adhesive systems as well as new 
types of brackets are being improved. Unfortunately, fixed orthodontic appliances cause a 
lot of doubts, because of their potentially destructive influence on the enamel [7,8] and 
others [9-12]. In [7], there are no statistically significant differences in volumetric changes 
after polishing between the different clean-up methods. However, sufficient clean-up with- 
out enamel loss was difficult to achieve. Enamel breakouts after debonding were detectable 
in 27% of all cases (mean depth of 44.9 mm). In [8], uncoated and precoated brackets exhib- 
ited similar debonding patterns. The debonding method tested in this study did not restore 
the original enamel surface, although there was no clinically relevant enamel damage. The 
use of coherence tomography in the areas of the mentioned enamel thickness or in the as- 
sessment of the strength of bonding is presented both in [9], with the use of a dental OCT 
handpiece, and in [10,12]. Therefore, there is a need for quantitative and not only qualitative 
assessment of tooth enamel in order to select the appropriate methodology and course of 
treatment as well as to assess the quality of enamel after treatment [13-16]. 

Evaluation of enamel quality can be done in many different ways. The known methods 
use an atomic force microscope (AFM) and a scanning electron microscope (SEM). It can 
also be evaluated using optical coherence tomography (OCT) [12]. This imaging method 
enables non-invasive structure visualization on the maximum enamel thickness and the 
desired tooth surface. The use of an optical tomograph having a wavelength of 840 nm to 
evaluate the enamel thickness is also interesting. For example, in [9] a central wavelength 
of 1310 nm and a spectral bandwidth of 47 nm were used for this purpose. 

Unfortunately, at present it is relatively difficult to compare images taken at different 
time intervals and at different stages of treatment (diagnosis) [17-19]. The difficulty lies 
in automatic finding of the same areas or points on two or more images. Therefore, 
there is a need to propose a new automatic algorithm able to analyse OCT images ac- 
quired at different time intervals and quantitatively assess the treatment progress 
[20,21]. A profiled automatic algorithm for the analysis and processing of OCT images 
of tooth enamel, in particular, should enable to: 

• quantify the changes in the enamel structure after treatment with fixed braces 
depending on the used orthodontic brackets, methods of attachment and tissue 
development after bracket removal, 
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• quantify the efficiency of cleaning a tooth surface off the composite material 
residue after the removal of orthodontic brackets when using selected ways of 
enamel development. 

A description of the dedicated automatic algorithm is discussed in detail later in 
the paper. 

Material 

OCT images of 180 teeth were obtained from the Topcon 3D OCT-2000 camera - Figure 1. 
The images were obtained in vitro by performing sequentially 7 stages of treatment on each 
of the 180 teeth: before any interference into enamel, polishing with orthodontic paste, etch- 
ing and application of a bonding system, orthodontic bracket bonding, orthodontic bracket 
removal, cleaning off adhesive residue. 7 OCT images were obtained at each of the 7 stages 
(the first image refers to a tooth without any interference). 

All the teeth were extracted and then stored for 30 days in distilled water with 0.1% 
crystal thymol at room temperature, and then cleaned with pumice using a eraser (Top- 
Dental), rinsed with distilled water and dried with compressed air for 15 seconds. The 
composite material BluGloo (ORMCO, USA) was used to attach orthodontic brackets. 
The vestibular tooth surface was etched for 30 seconds with 36% phosphoric acid - Blue- 
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Etch (CERKAMED, Poland), rinsed with distilled water and dried with compressed air. 
The adhesive system Ortho Solo (ORMCO, USA) was applied to the etched surface. A 
bracket was attached with the use of clamping tweezers for tooth enamel so that its centre 
was 3.5 mm below the edge of the occlusal surface, in the middle of the mesial-distal axis 
of the tooth. The test teeth were stored in water at 37 °C for 24 hours. Adhesive material 
residue was removed from enamel in two ways: by the use of a carbide cutter and by 
means of cup-shaped elastics OneGloss 0181 (Shofu Dental Corp., USA). 

The acquired images (B-scans) had a resolution of MxN= 884x512 pixels (where 
M - number of rows, N - number of columns) where a single pixel covered the area of 
5 x 11.7 (im at the colour resolution of 8 bits/pixel. For each patient /= 128 B-scans 
(every 47 |im) were obtained, which allowed for full 3D tooth image reconstruction. All data 
were analysed in a source format *.fda, *.fds. The data were anonymised and tooth extrac- 
tion in all the patients was performed in accordance with the Declaration of Helsinki. 



Method 

Pre-processing 

Pre-processing of images is related to automatic reading of a sequence of OCT images 
from a source file with the extension *.fds, thus obtaining a matrix of individual images. 
The read data also include information on the location of subsequent B-scans, their size, 
date of test, patient data and others. The sequence of images L GRA ^m,n,i) (Figure 2) thus 
obtained is further subjected to median filtering with a mask h 1 sized M hl xN hl x I hl = 
3x3x3. The mask size was chosen arbitrarily taking into account the size of noise and 
artefacts entering the optical path. The filtered image L M (yn,n,i) (where m,n are the coor- 
dinates of the rows and columns of the matrix me(l,M) and «e(l,A/)) is further subjected 
to preliminary transformations. These include determination of the enamel layer, meas- 
urement of its thickness in the entire area, comparison with other images of the same 
tooth and determination of the value of the proposed enamel quality coefficient. 



LGRAY(m,n,i) 



n [pixel] 



ax 



m [pixel} 



u <u <tt m 



n [pixel] 

Figure 2 Reconstruction of a sequence of images L GRAY {m,n f i) for M x N x / = 884 x 51 2 x 1 28 pixels. 

The image shows the result of reconstruction of an image sequence with a sample B-scan obtained for 
i = 60. The analysis aims at automatic determination of the edge on a sequence of OCT images (B-scans) in 
order to determine enamel thickness. 
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Processing 

The filtered image L M (yn,n,i) free of noise and artefacts is the basis for further transfor- 
mations. These include determination of enamel thickness and imposition of two 2D 
images (C-scans) in order to compare changes in enamel thickness. 

Determination of enamel thickness 

Determination of enamel thickness is related to determining the position of two enamel 
layers: outer L en (n,i) and inner L in {n,i) directly in the image L M (m,n,i). Determination of 
the outer layer L en {n,i) started from setting the upper limit of L e ^n,i) between the object 
(objects) and background: 

Lef(n, i) = mm(Lefrn(m, n, i)) (1) 

m 

where: 

L efm {n,i)-^ N Qther (2) 

and p r - binarization threshold determined automatically - Otsu [22]. 

The upper limit of the object L e j(n,i) is shown in green in Figure 3a. The presented 
method according to the equations (1) and (2) gives correct results in 100% of cases 
(for more than 20,000 images). This limit is used to determine the correct upper and 
lower boundaries of enamel (only a portion of the upper limit of objects L e j(n,i) is the 
appropriate outer boundary of the enamel L en (n,i)). 

Detection of the enamel lower boundary L in (n,i) is related to the analysis of each col- 
umn of the image L M {m,n,i) in the range from L e j(n,i) to the last row. Sample data 
(levels of brightness) are shown in Figure 3b and consist of two parts. Part I (Figure 3b 




n [pixel] 



Figure 3 Example of the result of the upper boundary analysis and selected parts of the analysis. 

The image a) shows the result (in green) of automatic detection of the upper boundary L ef (n,i). The image 
b) shows changes in brightness for the selected A-scan {L M {m,n = 240, /' = 60)) and the image c) shows the 
histogram on the basis of which the binarization threshold is selected in order to detect the enamel 
inner boundary. 

\ J 
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concerns reduction of brightness of a fixed speed whereas part II refers to a constant 
brightness level Part III is a shift of the mean (line II) by a standard deviation of the 
mean from this range to the top. The grey level in the range from p r2 to p r3 (II to III) is 
a place indicating tooth translucency reduction related to the end of the enamel layer 
and the start of another layer, namely dentine. The decision which of the threshold 
values should be taken as the correct one (in the range from p r2 to p r3 ) is connected in 
this case with a limit shift of about 70 pixels which corresponds to 70-5 = 350 um. The 
value of this threshold can be automatically determined and clarified on the basis of 
the previous and subsequent columns and histogram analysis - Figure 3c. The thresh- 
old value Pr2( n >i) (simplified to p r2 ) was calculated automatically based on the 
histogram: 



p r2 (n, i) = arg^max(L to («, b, i)) ) (3) 
where: 

M 

L hist (n, b, i) = ^L s (m, n, b, i) (4) 

m=l 

L s (m,n,b,i) = [ Q Qther (5) 

The obtained result, image L B (m,n,i) for the adopted binarization threshold, is shown 
in Figure 4a: 



t ( »\ / 1 if L M (m,n,i 

**(*.».*) = { 0 other 



)>Pr2 (6) 



The binary image L B (m,n,i) is further subjected to median filtering with a mask h 2 
sized Mh 2 x N^ 2 x 4 2 = 3 x 3 x 3. Then, the resulting image L MB (m,n,i) undergoes a mor- 
phological closing operation, i.e.: 

L c (m,n,i) = mini max(I M5 (m, n, i)) I (7) 

SE V SE J 

where SE - structural element sized 3x3. 

The size of the structural element was chosen on the basis of measurements and 
average noise sizes (isolated ones) that make up clusters consisting of no more than 
2 or 3 pieces. The next step involves filling holes and labelling the image L c {m,n,i). 
All further operations are carried out for each of the two-dimensional images for 
constant z, that is L c (rn,n) (L c (yn,n,i = const)). Labelling enables to give a separate 
label (number) to each cluster of ones L^m^). In this case, it is possible to elimin- 
ate all of the clusters which do not have a maximum surface area. It results in the 
image Lw{m,n) shown in Figure 4. In the next step, the inner enamel boundary L in {n) is 
determined as a 3rd order polynomial approximation of the course of the inner 
edge L inr {n) of an object in the image L^m^). An example of a polynomial takes 
the following form: 

L in {n) = -0.0001 -n 3 + 0.018-n 2 -3.21-n + 693.39 (8) 
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b) 




E 




Figure 4 Example of the binarization result of an OCT image and the results of analysis of inner 
and outer boundaries of the enamel. The image a) shows problems that occur when determining the 
enamel inner boundary on a binary image whereas the image b) shows an example of the result of 
analysis based on the proposed algorithm. The outer boundary is marked in green, the inner boundary 
before filtration in red and the inner boundary after filtration in blue. 



The outer edge of the enamel, L en (n) for short, is part of the previously defined edge 
of the object L e fn) (for i = const), i.e.: 

L (n) = i Le ^ n) if (Lef(n) + 1 = 1) A(L ef (n)<N) {g) 
\ 0 other 

Based on the determined edges L e j(n) and L in {n), their difference is calculated: 

L d {n) = L in (n)-L en (n) (10) 

This difference is calculated for every i (for each B-scan), i.e. L d (n,i) forms a matrix of 
differences. The obtained matrices of differences for individual tooth treatment steps are 
shown in Figure 5 in the following order: the tooth before interference L da (n,i) - Figure 5a, 
after polishing with paste L db (n,i) - Figure 5b, after etching L dc (n,i) - Figure 5c, after the 
application of the adhesive system L dd (n,i) - Figure 5d, the orthodontic bracket L de (n,i) - 
Figure 5e, after orthodontic bracket removal L d pi,i) - Figure 5f, after cleaning off adhesive 
residue L dg (n,i) - Figure 5g. These matrices (from L da (n,i) to L dg (n,i)) are considered in the 
subsequent stages of analysis. 
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n [pixel} 



Figure 5 Examples of results of enamel thickness measurement for the subsequent stages of tooth 
treatment. It shows: a) enamel thickness before any interference L da {n,i), b) after polishing with paste L db {n,i), 
c) after etching L dc {n,i), d) after application of the adhesive system L dd {n,f), e) an orthodontic bracket L de (n,i), 
f) after bracket removal L d in,i), g) after cleaning off adhesive residue L dg {n,i). These images then form the basis 
for the analysis of the impact of treatment methods and choice of orthodontic systems on changes in enamel 
thickness. The thickness is expressed in pixels for the artificial colour palette (ranges from blue to red, and 
passes through the colours cyan, yellow, and orange). 
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Imposition of two 2D images 

The images from L da (n,i) to L dg (n,i) show enamel thickness in various stages of tooth 
treatment. From a practical point of view, a change in thickness is important and, to be 
more exact, a change in enamel volume for each step. Therefore, it is necessary to 
propose a method of automatic imposition of individual images from L da (n,i) to L dg (n,i). 
Imposition of images, finding homogeneous points, is connected in this case with finding 
three values, namely a shift in the axis An, a shift in the axis Ai and zooming Az. In prac- 
tice, due to the specificity of attachment used, the errors resulting from rotation did not 
exceed 1°. Because of this, it was assumed that there is no error in the images caused by 
rotation. The searched shifts in the axes An and Ai are further expressed in pixels whereas 
zooming Az is expressed as a percentage of the original image, i.e. 10% means a 10% zoom 
with respect to the reference image. 

It is most convenient to find homogeneous points in standardized input images after 
median filtering and not in images showing different enamel thickness, i.e. (assuming 
that the denominator is different from zero): 



The first image L Ta (n,i) is here treated as the reference one, symbols of the other ones 
are analogous to L d (n,i), these are subscripts from l V to ^ i.e. from L Tb (n,i) to L Tg (n,i). 
For each new shift and zoom relative to the reference image L Ta (n,i), the value of the 
match criterion / was calculated. The criterion / is the sum of absolute differences be- 
tween the reference image L Ta (n,i) and the image being compared, e.g. L Tb (n,i). There- 
fore a point reaching the minimum and the values of An, Ai and Az for which it 
occurred are looked for. It is expressed as a percentage of changes with respect to the 
reference image, i.e., e.g. -10% means 90% of the reference image for a given magnifica- 
tion of the image (superscript is expressed as a percentage, e.g. for 120% relative to a 
standard, L Tb 100 \n + An,i + Ai). Therefore, the criterion / calculated for a pair of images 
L Ta and L Tb takes the following form: 



for A«e(-40,40) pixels and A/e(-40,40) pixels. 

The missing or redundant rows and columns in relation to L Ta (n,i) were removed or 
zero was written so that the numbers of rows and columns are identical. The position 
in both axes was changed every pixel, whereas zooming was carried out in the range 
from -20% to 20% every 10% relative to the reference image. The zoom for which cal- 
culations were made was marked as a superscript, e.g. if there is no image magnifica- 
tion then L Tb °\n + An,i + Ai). The obtained results, changes in the criterion / for 
different positions of the image L Tb °\n + An,i + Ai), in the absence of magnification, 
are shown Figure 6. After imposition of all the images and appropriate adjustment of 
the position of images from L db (n,i) to L dg (n,i) (after adjustment starting from L db \n,i) 
to L dg (n,i)), a region of interest (ROI) is determined. It has to meet two criteria: its 
total area must be within the matrix area and it must coincide with the contour of the 
tooth. The second criterion requires determination of the tooth area. 




(ii) 




i=l n=l 



(12) 
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Figure 6 Graph of changes in the criterion J for a sample pair of images L Tb (0) {n + An,/ + A/) and L Ta {n,i) 
and images after correction. The value looked for in the graph is the minimum for which the values of shifts 
in both axes (An, A/') and the zoom which gives the smallest value are read. In this case, the graph refers to a 
constant zoom of the image L Tb (0) (0%) and for the example shown the following values were read from the 
graph: An = 9 and A/' = 4 pixels. In the left hand side of the image, there are the results obtained after shifting: 
images L Ta and L Tb *. 



Determination of the tooth contour 

Determination of the tooth contour is related to the analysis of one selected image L da 
(n,i) or L Ta (n,i). The other images from L db \n,i) to L dg \n,i) have contours in the same 
place, which results from their earlier imposition (finding homogeneous points). Just 
like in the case of imposition of images, the image L Ta (n,i) is easier to analyse than its 
corresponding image of the enamel thickness L da (n,i). In the first step, the image L Ta (n,i) 
is subjected to median filtering with a mask h 3 sized N^ 3 x I h3 = 31 x 31 pixels. The mask 
size was selected on the basis of the average width of the tooth and the size of artefacts 
being the image noise. The output image L Ra (n,i) is the basis for determining the contour 
using the standard method and the active contour method. For the latter method, the fol- 
lowing operating parameters were adopted: changes in the axes of rows and columns in 
the range of ±20 pixels and averaging the grey level in the areas sized 6x6 pixels on both 
sides of the analysed pixel. The algorithm stops work at a time when there is no evidence 
of differences between successive steps of bringing the contour closer to the object. The 
results are shown in Figure 7. The determined ROI is the basis for further analysis of the 
results. The resulting images from L ROIa to L ROIg contain information about the thickness 
of the layer only in the designated area of the tooth. 
A block diagram of the proposed algorithm is shown in Figure 8. 



Results 

The obtained images from L ROIa to L ROIg containing information about the thickness of 
the enamel layer are the basis for the calculation of mean thickness, standard deviation, 
and minimum and maximum values. The results are shown in Table 1, whereas Table 2 
presents information about enamel thickness differences between the images. According 
to the adopted symbols from L ROIa to L ROIg , these are: tooth before interference 
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Start Z-Ra 

\~End\ 




Figure 7 Figure showing operation of the active contour method. The input image L Ta after median 
filtering (L Ra ) is analysed using the active contour method. The standard whose size is independent of the 
size of the object is highlighted in green in the image. The blue shows the final stage of approximation - the 
area of ROI analysis. 



(index 'a'), after polishing with paste (index ( b'), after etching (index ( c), after application 
of an adhesive system (index l d£\ an orthodontic bracket (index V), after removal of the 
orthodontic bracket (index '/) and after cleaning off adhesive residue (index The dif- 
ferences confirm enamel thickness loss of 80 \im (enamel thickness at this stage is 650 ± 
129 (im) after polishing with paste, loss of 435 \im (enamel thickness at this stage is 295 ± 
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Figure 8 Block diagram of the proposed method for image analysis and processing. First, the input 
image is subjected to pre-processing (reading of data from the format *fds, filtration) and then the enamel 
layer is recognized on subsequent B-scans. After that, C-scans obtained for subsequent phases of the enamel 
treatment are imposed on each other and the difference in enamel thickness in the ROI is measured. The last 
stage of processing involves comparison of changes in the thickness of enamel - its volume for the subsequent 
processing steps. 
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Table 1 Mean, minimum and maximum values of the enamel thickness in the ROI and 
standard deviation of the mean for the subsequent processing steps (according to the 
ROI designated in various stages - Figure 5) 

Subsequent processing steps [\xm] 



Parameter 


l-ROIa 


l-ROIb 


Lroic 


LROId 


l-ROle 


l-ROIf 


l-ROIg 


Mean 


730 


650 


295 


560 


80 


455 


605 


std 


165 


129 


55 


98 


135 


91 


140 


Min 


390 


355 


125 


25 


5 


5 


240 


Max 


1105 


1015 


515 


820 


775 


665 


860 



55 (im) after bonding, an increase in the layer thickness of 265 um (enamel thickness at 
this stage is 560 ± 98 um) which is the bonding system. After the removal of the ortho- 
dontic bracket, the adhesive residue is 105 um (enamel thickness at this stage is 455 ± 
91 um) and after removing the adhesive residue, the enamel thickness is 605 um. The 
thickness of the enamel before and after the whole treatment decreased by about 125 um. 

The results in Table 1 contain information about the minimum and maximum tooth 
enamel thickness. The differences in these values (minimum and maximum) result 
from actual variations in the enamel thickness in various locations of the tooth and the 
method errors. These errors are influenced to the greatest extent by the automatic se- 
lection of binarization thresholds p r2 and p r3 shown in detail in Figure 3 and the conse- 
quences in the binary image in Figure 4. 

Accuracy verification of the proposed method 

The accuracy of the proposed method was assessed on the basis of the comparison of the 
results of enamel thickness measurements obtained with the proposed algorithm with in- 
dependent measurements carried out by an expert. They involved manual marking of the 
enamel boundary by an expert on a tomographic image for 20 teeth (20-128 = 256 B- 
scans). In the case of the manual method, an expert dentist indicated two boundaries: the 
enamel boundary L in (E \n) and the outer tooth boundary L e J E \n) - Figure 9. On this basis, 
the enamel thickness was determined (analogously to (10)): 

Lf(n)=L^{n)-L^{n) (13) 



Table 2 Changes in the enamel thickness in the ROI for the subsequent processing steps 
(according to the ROI designated in various stages - Figure 5) 
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Figure 9 Results of the automatic operation of the algorithm and the manual method of determining 
the enamel thickness. The expert dentist marked two contours - the outer tooth contour L en (E) (n) (in magenta) 
and the enamel boundary L in (E) (n) (in yellow). On this basis and on the basis of the results of automatic analysis 
{Lenin) green, L in {n) blue) with the proposed algorithm, the measurement error was calculated per 20 teeth. This 
error did not exceed 15%. 



The measurement error, adopted as the measurement error of the proposed algo- 
rithm, was designated as: 



1 » \U(«)-Lf(n) 



For the analysed 20 teeth, the error 8 d did not exceed 15%. Individual results of the 
manual and automatic analyses are shown in Figure 9. In none of the analysed 20 teeth, 
the algorithm was mistaken by more than a few pixels when defining the outer tooth 
layer - L en (n). The contour of this layer was in each case very clearly visible and the dif- 
ferences between layer designation by an expert and the algorithm were visible only 
after zooming (in the lower part of the image Figure 9). In contrast, the enamel layer 
marked by the expert was only visible in selected ranges of columns. This situation, for 
layers L in (ri) and L in (E) {ri), is shown enlarged in Figure 9 at the bottom. Individual yel- 
low points spaced apart by a few pixels indicated by the expert are visible. These points 
were further connected with segments. It was possible owing to information derived 
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from anthropometric data, on whose basis the curvature of the contour and its possible 
variation per a single tooth are known. In summary, using the described algorithm, the 
measurement error b d does not exceed 15%. 

Comparison with other known methods 

The method of analysis and processing of OCT images of tooth enamel developed by 
the authors can be compared to other known methods of analysis. These include: 

• methods of analysis and processing of OCT images applicable in ophthalmology 
(e.g. recognition of retinal layers) and 

• teeth image analysis carried out in visible light. 

In the first group, there are methods which enable detection of individual layers in 
OCT images of the eye fundus [23-25]. The class of images is very similar to OCT im- 
ages of tooth enamel. These methods work in a time comparable to that of the algo- 
rithm (from 10 to 200 ms/2D image) [26-28]. However, specificity of OCT images for 
the eye enables a significant simplification of the algorithm in the detection area of in- 
dividual layers. This is due to a much better contrast between the layers. The advan- 
tage of OCT images of tooth enamel is the necessity of determining usually only one 
layer. 

The other group of methods is closely related to dentistry but the ray spectrum range 
is significantly different (visible light). The acquired images do not allow for assessment 
of the thickness of any tissue, however, some of the methods of image imposition can 
be transferred to OCT imaging. An important difference is the possibility of using, in 
the case of visible light, RGB components which enable to distinguish a tooth from a 
handle and background. The applied active contour method brings correct results in 
both cases (visible light and OCT). 

The forms of imaging, namely in vitro and in vivo, are also extremely important here. 
The differences are very interesting not only from the clinical point of view but also in 
terms of the algorithm of image analysis and processing. In the first case (in vitro), re- 
peatable attachment of the tooth greatly facilitates image analysis - the segmentation of 
the tooth boundaries. In the other case (in vivo), the situation is more difficult because 
not only the background but also the mutual adjacent teeth arrangement change. In 
this case, there may also occur other problems with imaging such as placing the tomo- 
graph perfectly perpendicular to the tooth etc. However, preliminary tests confirm the 
usefulness of the proposed method, also in in vivo studies. The proposed segmentation 
process is resistant, in the area declared, to the impact of the background of the sur- 
rounding teeth (in the case of in vivo tests) on the results obtained. This is achieved 
owing to a segmentation method using the active contour. In this method, an import- 
ant element is the contour of the tooth and not the background as may be the case in 
the method based on thresholding (binarization). 

The presented method can also be compared to other methods of imaging of the en- 
amel layer. As mentioned in the introduction, the known methods of tooth enamel ana- 
lysis include assessment by means of an atomic force microscope (AFM) [29,30] and a 
scanning electron microscope (SEM) [31]. The class of images and the depth of im- 
aging in these cases are quite different from the ones in OCT and typically confined to 
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the outer enamel surface analysis - like, e.g., in [29]. A different class of images deter- 
mines a different type of a profiled algorithm for image analysis and processing. 

The other known methods of enamel thickness analysis do not enable automatic, 
quantitative measurement of the enamel thickness present in the ROI and automatic 
comparison of image groups. This is the case in [32], where comparisons between 
specific areas of the tooth enamel were made manually in OCT images. Automatic 
measurement was presented only in [33]. However, it concerns polarization sensitive 
optical coherence tomography (PS -OCT) and is not related to the problem of overlap 
of individual images in the subsequent processing stages of the tooth as shown in this 
paper. 

Summary 

The paper presents the possibility of using the profiled algorithm proposed by the au- 
thors for the analysis of OCT images of tooth enamel. The algorithm of OCT image 
analysis and processing has the following features: 

• works fully automatically - without any operator intervention, 

• allows to quantify the changes in the structure of enamel, 

• allows for quantitative assessment of the effectiveness of cleaning the tooth surface 
and the effectiveness of the use of selected methods of enamel development, 

• analysis time of a sequence of 2D images does not exceed 5 seconds when using the 
Core i5 CPU M460 @ 2.5GHz 4GB RAM, 

• the results of the mean thickness of the tooth enamel, minimum and maximum 
values as well as standard deviation are analysed automatically and saved to text 
files *.txt and Excel *.xls. 

The presented algorithm is only one possibility of completing the task described in 
the paper. Automatic analysis of tooth enamel thickness provides a number of further 
possibilities. These include area analysis of enamel thickness (for each individual tooth 
area separately) and enamel texture analysis. Imaging and quantitative measurement of 
the enamel structure before installation of braces and after their removal enables to ex- 
pose the tooth tissue damage extent depending on the used brackets and method of at- 
tachment. This makes it possible to deduce which brackets and what technique of their 
installation is the safest for tooth enamel. 

The presented methodology does not fully exhaust the issue of automatic measure- 
ment of enamel thickness and the effect of the method of attachment of orthodontic 
brackets. The presented tool will be used in future work to assess the impact of mul- 
tiple parameters. For example, the impact of the OCT settings on the results will be 
analysed. Moreover, the effect of time of in vitro studies on the result will be assessed. 
The analysis of reproducibility of results obtained in tomographic cameras at different 
wavelengths in the range from 600 to 1400 nm is equally interesting. An important as- 
pect here will be confirmation of the universality of the proposed algorithm of image 
analysis and processing for different types of OCT devices or even the impact of the 
type of studies, in vitro or in vivo, on the process of segmentation. 

It is also possible to modify the approach to the presented problem by using, for ex- 
ample, texture analysis [34], morphological approach [35], spectral analysis or other 
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methods used for slightly different applications which, in practice, generate correct results 
[36-41]. In the future, the authors intend to perform and describe these methods, espe- 
cially area and morphological analysis of the tooth enamel, by comparing the results ob- 
tained and measurement errors resulting from various adopted methods of analysis. 
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